UNCLASSIFIED 


ad  429850 


DEFENSE  DOCUMENTATION  CENTER 

FOR 

SCIENTIFIC  AND  TECHNICAL  INFORMATION 

CAMERON  STATION.  ALEXANDRIA.  VIRGINIA 


UNCLASSIFIED 


NOTICE:  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  incurs  no  responsibility,  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 


COLLAPSE,  BUCKLING  AND  POST  FAILURE 
OF  CYLINDRICAL  SHELLS 


o 

Q 

1 

Q 

t 

>* 

J 

o 

CD 

LlJ 

CD 

CD 

i _ _ 

<£ 

f 

AS 

Joshua  E.  Greenspon,  Dr.  Eng. 


O 

10 

00 

O) 

CQ 


Contract  No.  DA  36-034-ORD-3081  RD 
Technical  Report  No.  4 
December,  1963 


-Tr 


J  G  Engineering  Research  Associates 


3831  MENLO  DRIVE 


BALTIMORE  15,  MARYLAND 


COLLAPSE,  BUCKLING  AND  POST  FAILURE 


OF  CYLINDRICAL  SHELLS 


Joshua  E.  Greenspon,  Dr.  Eng. 

J  G  ENGINEERING  RESEARCH  ASSOCIATES 
3831  Menlo  Drive 
Baltimore,  Maryland  21215 


Ballistic  Research  Laboratories 
Aberdeen  Proving  Ground 
Contract  No.  DA  36-034-ORD-308I  RD 
Technical  Report  No.  4 
December,  1963 


(Reproduction  in  whole  or  in  part  is  permitted  for 
any  purpose  of  the  United  States  Government) 


ABSTRACT 


In  this  report  it  is  demonstrated  how  simple  existing  theories  can  pre¬ 
dict  whether  a  shell  will  fail  by  buckling  or  form  a  collapse  hinge  and 
then  fail  by  this  collapse  mechanism.  It  is  also  shown  how  work-energy 
analysis  can  predict  the  post  failure  plastic  deflections  of  the  shell 
once  the  mode  of  failure  is  known.  Some  comparisons  with  experimental 
results  are  given. 
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stress  resultants 
longitudinal  coordinate 
peripheral  coordinate 

arbitrary  functions  of  (ft  which  are  determined  from  the 
boundary  conditions 

mean  radius  of  the  shell 

pressure  distribution 

thickness  of  shell 

longitudinal,  tangential  and  radial  displacements  respectively 
parameters  used  to  define  the  peripheral  pressure  amplitude 
Poisson's  ratio 
modulus  of  elasticity 


JL 
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CfX 

(%Cr)b 

f 

V 

sur6 

<Zc 


length  of  shell 

•umber  of  circuit: 'ierential  waves  in  buckling  pattern 
yield  stress  in  pure  tension 
pressure  which  initiates  collapse 
pressure  which  initiates  buckling 
mass  density  of  shell  material 
spatial  distribution  of  the  impulse 
work  done  by  shell  in  deforming 
maximum  deflection  amplitude 
spatial  distribution  of  the  deflection 
yield  strain  in  pure  tension 
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K 

TJ 

x 

£e 

A 


*/JL 

width  of  plastic  hinge 

defined  by  Eq.  [20] 

total  impulse  (theoretical) 

defined  by  Eq„  [22] 

experimental  impulse  per  unit  area 

reflected  experimental  impulse  per  unit  area 

incident  experimental  impulse  per  unit  area 

parameter  defining  periphery  die  out  of  spatial  deflection  dis¬ 
tribution  in  buckling 


I.  Introduction 

1* 

A  previous  report  presented  the  plastic  work  expressions  for  cylindri¬ 
cal  shells  and  discussed  the  types  of  failure  that  could  o^cur  in  cylin¬ 
drical  shells  under  blast  loading.  Experiments  have  shown  that  cylin¬ 
drical  shells  subjected  to  side-on  blast  can  go  into  two  main  types  of 
failure.  These  are  buckling  and  collapse.  The  buckling  type  of  fail¬ 
ure  is  described  by  a  deformation  pattern  which  consists  of  a  number 
of  lobes  around  the  periphery  of  the  shell  and  one  half  wave  length 
along  the  length  as  shown  in  Figure  la.  The  collapse  failure  is  de¬ 
scribed  by  a  straight  failure  hinge  at  the  center  of  the  shell  as 
shown  in  Figure  lb.  ^oth  of  these  figures  are  taken  from  Schuman's 
experimental  results.  The  type  of  failure  will  depend  upon  the  geo¬ 
metry  of  the  shell  and  can  be  predicted  from  an  elastic  stress  and 
buckling  analysis  of  the  shell  as  will  be  seen  later  in  this  report. 

II  Elastic  stress  analysis 

Assume  that  the  shell  is  thin  and  that  membrante  theory  is  adequate 
to  describe  the  stress  patterns  in  the  shell.  Assume  also  that  the 
shell  is  of  length  and  is  supported  at  each  end.  Take  the  origin  of 
coordinates  at  the  center  of  the  shell  as  shown  in  Figure  2.  The  mem¬ 
brane  forces  are  shown  on  the  differential  element  in  Figure  3. 

If  •/pCz-jCp)  is  the  static  load  per  unit  area  applied  la^e^ally  to  the 
shell  then  it  follows  from  the  basic  membrane  equations  '  that 


Ny  -  (X  ZjCf) 


where  t  <V)  and  are  functions  of  <P  which  are  to  be  determined 

from  the  boundary  conditions  on  Nxcf  N cp  '  if  some  °f  the 

boundary  conditions  are  given  in  terms  of  displacements  then  the  follow¬ 
ing  membrane  equations  in  terms  of  displacements  must  be  utilized 


*Superscripts  refer  to  references  given  at  the  back  of  the  report. 
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[2] 


EJlu  ~f(Njc  -i)hJc/)dx.  -h^Ccf) 

eAv  =  2.(i+J)JrN*.epc/x.  -  j-zf+{c/>) 

eA^=  fi-CNf-M*)*  EA^ 

where  A  i-s  the  shell  thickness,  £  is  the  modulus  of  elasticity, 
are  the  displacements  (see  Fig.  3)  and  E^Ccf) 

are  arbitrary  functions  to  be  determined  from  the  boundary  conditions. 

Now  assume  that  the  pressure  'PfejCf)')  can  be  represented  as 

-  /£>  C  ^ o<ef*.7X 

-ftflK-cp)  [3] 

r  c  Ci7T-cP )  7C<cf<zrr 


where  -A  is  the  maximum  amplitude  of  the  pressure  and  /2> J  ^  are 
parameters  determining  the  peripheral  pressure  distribution.  If  the 
shell  length  is  small  compared  to  the  distance  from  the  explosion  and 
if  the  blast  is  a  side-on  load  then  it  seems  reasonable  that  the  spatial 
load  distribution  given  by  [3]  should  resemble  the  spatial  distribution 
produced  by  the  blast. 


The  solution  will  be  obtained  for  the  boundary  conditions 

u  =  o  «+  £  =  -  &/z. 


[4] 


By  straight  forward  integration  of  equations  [1]  and  [2]  subject  to 
boundary  conditions  [4]  it  is  found  that 


A/^>  -  (X  j~Ccf)  j  -rfy))-  C  [5] 

N,t--p^  ft 

N^-fs.  ^Pfcf)  +-)a.  ?.£<?)-  A 

The  stress  resultants  at  ZL  =  =  CD  are 

Ntf 

My  =  O  l6' 

Nx.  =  Oa -pc  - 

Yielding  will  occur  at  X.  -  cf  -  c>  when  the  Von  Mises  Yield 

Conditon  is  satisfied  at  this  point,  i.e. 

[7] 
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where  (77  is  the  yield  stress  in  pure  tension.  Substituting  the 
values  of  the  stress  resultants  into  the  yield  condition  it  is  found 
that  yielding  will  occur  at  the  center  of  the  shell  (  2L  -  o  ^  cfl  -  O  ) 

when  the  pressure  ^pa  has  a  value 


(■fX  -  °°  £  7  - - - 18] 

III.  Elastic  buckling  analysis 

The  classical  theory  of  buckling  of  cylindrical  shells  is  presented 
by  Timoshenko , 5  It  is  found  that  the  lateral  pressure  at  which  buck¬ 
ling  under  uniform  loading  will  occur  is  given  by  the  following  re¬ 
lation^ 


t  (yi  /  + 


2  -K 

/  + 


J  J 


19] 


where  7V  is  the  nurriber  of  full  waves  around  the  periphery  of  the 
shell.  It  has  been  found  by  Reynolds  that  in  lobar  buckling  such 
as  this  the  circumferential  parameter  ru  approximately  satisfies  the 

relation 


1 


I  -t- 


<JZi 


(10) 


w-  ^  7'  2  2  '"X 

and  that  buckling  will  always  occur  with  only  one  half  wave  along 
the  lenguh.  In  other  words,  if  Yn,  denotes  the  number  of  axial  half 
waves  along  the  length,  then  'yn.  =•  /  and  is  determined 

from  equation  [10].  The  factor  1.23  has  to  be  adjusted  so  that 
turns  out  as  a  whole  number.  Once  'Yi*  is  determined  then  equation 
[91  will  give  the  value  of  the  buckling  load  for  uniform  loading.. 

For  nonuniform  loading  correction,  factors  can  be  applied  in  accor¬ 
dance  with  a  recent  paper 


IV.  Buckling  or  yield 

It  is  clear  from  Schuman's  experiments2  that  both  buckling  and  yield 
collapse  can  occur..  The  main  problem  is  to  be  able  to  predict  which 
type  will  take  place.  Once  buckling  or  collapse  has  commenced  the 
plastic  deformation  will  take  place  in  that  particular  pattern  into 
which  the  failure  has  started.  Iu  will  now  be  shown  how  equation  [8] 
and  [9]  can  be  employed  to  determi.ne  whether  a  collapse  or  buckling 
will  take  place.  Although  there  were  some  gross  assumptions  involved 
in  the  derivation  of  equation  [8]  the  main  item  to  be  recognized  is 
that  the  pressure  for  this  type  of  yield  is  proportional  to  /(X.  and 
dependent  on  0>  &/&_  .  On  the  other  hand  the  buckling  pressure  is 

critically  dependent  on  'Yl  ^  and  *A  *"/ 3  as  well  as  >7  /ck 

I.t  should  be  made  clear  at  this  point  that  equation  [8]  and  [9]  are 
not  to  be  used  to  approximate  any  of  the  dynamical  parameters  of  the 
shell-  they  are  to  be  used  only  to  determine  which  type  of  failure 
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will  start.  Once  this  is  known  then  the  plastic  analysis  as  given  in 
a  previous  reference^  will  be  employed  to  determine  the  plastic  deflec¬ 
tions  and  impulse  values.  Assuming  that 

£  X  /ooo  ^  j  /?=  2.^  °<  =  2 

it  was  found  that  shell  buckling  or  yield  could  be  predicted  by  the 
equations  in  the  previous  section  by  employing  the  following  criterion: 

If  the  buckling  load  is  less  than  the  yield  load  the  shell 
should  buckle;  if  the  yield  load  is  less  than  the  buckling 
then  the  shell  should  collapse. 

Due  to  lack  of  experimental  evidence  on  the  load  distribution  the 
parameters  and  are  questionable.  Perhaps  when  more  information 
is  available  about  the  load  distribution  more  accurate  computations 
can  be  performed.  However  our  present  calculations  do  indicate  that 
the  simplified  equations  presented  above  describe  the  physical  phenom¬ 
enon  and  enable  prediction  of  the  type  of  failure  that  will  occur  in 
a  given  cylindrical  shell. 

To  illustrate  the  procedure  for  determining  whether  collapse  or  buck¬ 
ling  failure  will  occur,  consider  several  cylindrical  shells  that  were 
tested. 

1.  Consider  a  steel  shell  with  radius  of  1.5", length  11.62",  thick¬ 
ness  .019"  subjected  to  a  lateral  blast  which  is  8'  away.  The 
experiment  showed  that  this  shell  formed  a  collapse  hinge. 

Using  an  E  =  1000  =  2 ,  =  h,  0  =  .3  the  collapse  load 

according  to  [8]  is 

(^Pc)c  ~  07 

If  the  shell  would  buckle  [10]  predicts  that  0^.-3,  Then  equation 
[9]  shows  that  the  uniform  buckling  load  for  this  'Y-  is 

(fcr)h  a  ,OC2.l  <To 

and  using  the  nonuniform  pressure  correction  it  is  seen  that  this 
value  of  .0021  Cl  will  be  somewhat  higher,  possible  as  much  as 
50%  greater.  Thus  for  this  shell  ('pt)c  so 

the  shell  can  be  expected  to  fail  by  collapse. 

2.  Consider  next  an  aluminum  shell  of  radius  1.5",  length  6",  thick¬ 
ness  .006"  subjected  to  a  lateral  blast  which  is  15'  away.  The 
experiment  showed  that  this  shell  buckled  with  6  lobes  around  the 
periphery  and  a  half  wave  along  the  length.  Using  an  E  =  1000  <7 

-  2,  ^  =.3,  the  yield  load  would  be 

fPo  1  Cl 

Under  buckling  [10]  predicts  'K  -  £  ,  Equation  [9]  shows  that 
the  uniform  buckling  pressure  will  be 

(f-c  r)b  ^  •OOOZ^f-  Cl 

The  correction  for  nonuniform  loading  could  possible  raise  this 
value  to  .00036  67  .  Thus  for  this  shell  (fycr)  b  <  (po)c.  so 

that  the  shell  can  be  expected  to  buckle . 
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V.  Post  failure  collapse  and  buckling 
A.  Collapse 

It  has  been  found'*'  that  the  relation  between  the  total  impulse  applied 
to  the  shell  and  the  work  done  on  the  shell  can  be  written  as 


X‘ '/. dA 


[11] 


where  ft’*-.*')  is  the  spatial  distribution  of  the  impulse,  f 

is  the  mass  density  of  the  shell  material,  A  is  the  shell  thickness, 
and  \/  is  the  work  done  on  the  shell  in  deforming  it  plastically. 

The  integral  is  taken  over  the  ajea,  A  ,  of  the  shell.  By  a  systematic 
simplification  it  has  been  shown  that  V  can  be  written  as 

\ /„  zcr0Aisj-ajL 


AC z'jCP)cJ4>c/-)l 


\J  3  S0  Sq 

where  C To  is  the  yield  stress  in  pure  tension,  4^Ta  is  the  maxi¬ 
mum  deflection,  A  is  the  thickness,  A  is  the  length  and  TCx'ef) 
is  the  spatial  distribution  of  the  deflection,  i.e. 


[12] 


sur-  a*tc  ffiCjCp) 


The  assumptions  under  which  [12]  was  derived  are 

1.  The  shell  is  made  of  a  perfectly  plastic  material 

2.  The  radius  of  the  shell  is  considerably  smaller  than  the 
length  ( <</J 

3..  The  deflection  of  the  shell  is  much  smaller  than  the  radius 
(  ***/*.  <  <  /  ) 

where  is  the  yield  strain 

For  details  see  Reference  1 

The  deflection  pattern  for  collapse  is  shown  in  Figure  4  and  can  be 
written  analytically  as  follows 


AArC^cp)  -  a  dxxcf  -  \/q  a  -  if  Ldc  -  ^  xj  ^ 

Letting  x.1  -  Z-/4. 

W  ( '*.,(()  =  CX  Ceicfi  -  \/i3  *■-  siy,-  2 for  0<X.'<  -L 

=  a  ctfrcp- 


[13] 


[14] 


<0 


where  Jo  is  the  width  of  the  hinge  line  as  shown  in  Figure  4.  Using 
these  deformation  expressions  and  [12],  the  work  done  on  the  shell  in 
deforming  it  plastically  can  be  written  as 
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where 


2.*. 


<^-  7  C'+lX.')  fa^  -l  <JC  <•  o 

<2*.  * 

After  some  mathematical  manipulation  and  substitution  of 


it  is  found  that  \/  can  be  written 

'de/i 


\/  =  2Q-c&£  f,  , - -  ’  -/  \  / 

/J  ^  l  ^_v/  ?  ^ 

7  -  ^5^  / ] f,  _  )*■  7 

73  ds-LiKd)  4  z  5"/'  LtJ  J 


or  finally 


where 


Tf  1  — / 


<s/  “2^ 


[15] 


[16] 


[17] 

[18] 

[19] 


Now  using  the  impulse  equation  [11]  and  assuming  that 

The  relation  between  the  impulse  and  the  deformation  can  be  written 
as 


X,  ~  2jUci  Jk  1/  °JLl 
h  17  i/j 


[20] 


where 


i-.  _Z_  f  z  /dci-  M*"  cf  /  Z7TT5-  7 

*  =  ^/fC  TV- f -f  ?]/>-$)  J 


The  deflection  is  actually  described  by  do/d  .  A  plot  of 
as  a  function  of  d°/d  is  given  in  Fig.  5. 


Some  calculations  were  run  on  several  of  the  tested  shells  which 
failed  by  collapse  and  these  computations  with  the  experimental  re¬ 
sults  are  given  below.  A  given  value  of  deflection  corresponding 
to  a  given  do/d  was  measured  from  the  tests  and  the  incident  and 

reflected  impulse  per  unit  area  corresponding  to  this  do/d  were 
also  measured.  The  theory  (eq.  [20]  above)  predicts  thp  impulse  which 
will  produce  a  given  do/d  .  Unfortunately  the  exact  value  of  the 
experimental  impulse  was  not  given  but  the  incident  impulse*per  unit 

*The  incident  impulse  is  the  measured  impulse  without  an  obstruction 
such  as  the  shell.  The  reflected  impulse  is  the  impulse  that  reflects 
off  of  a  large  rigid  plane. 
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area  (T<)  £  ;  and  the  reflected  impulse  per  unit  area  ) r 

were  the  parameters  that  were  measured.  The  actual  impulse  applied__to  the 
shell  will  he  somewhere  between  these  values  of  and  (2Tt)r 

Assuming  that  the  measured  values  would  correspond  to  maximum  values  on 
the  shell  and  assuming  an  distribution  of  the  impulse  on 

the  shell  it  is  seen  that  the  total  impulse  would  be 

.  />  rr 


C^t)e=  if  f  C h)  £ 


(z)e 


~3Z 


So  the  theoretical  value  to  compare  with  the  measurements  is 

§=  /W 


A  value  of  C7J  -X 

shells  that  were  tested 


•SZJ  oc>c>  Ps  l 


was  assumed  for  the  steel 


Example  1 

3'f  Jz  ^  £  -  'Mi" 

Fcr-  de/cj'Zj  -frcr*  S~ 

S<s  -T^/23  psi  ~,l/is<c, 

The  measured  values  which  gave  the  shell  a  deformation  in  which  dt/d  ~  / 
(-Fe)l  -  psl  r*)lhs<£C,  j  fa),  -  29*9  p si  n*>l/iSe c. 


Example  2 


J  S  -  "  J  JLz,  C>9" 

Fc*  do/j  ~  ^  i//p  s i,  24 

X  cc.  2.  psi 


The  measured  values  were 


)i  ~  psi  r*  ilh see  .  ^  r  \JT 


Example  3 


d  -  3"  J  4*  r.a"  J  &  -,£>/<?" 

For  dc/gj  7>r  /7 z  1) 

-T  r  y/.sTpst 


The  measured  values  were 


(-F  a) i  &~j  (l-e.  }r  -  /  6>  O 

Example  4^  ^  £■=//,£  2."  >  h  -  ,  6,9  " 

her  do/j  7J-  ~  ,  2/  J“  ~  72.  ^~psC 

The  measured  values  were 

(rt);-:/+.r  /rt) 


»~il 4  sac 


-  SZ.  / 


The  different  impulse  values  arose  from  various  explosive  charge  weights 
and  distances  of  explosion  from  the  shell. 

The  comparisons  between  the  theoretical  predictions  and  experimental  re¬ 
sults  show  that  the  simplified  equations  predict  the  correct  order  of  magni¬ 
tude  for  the  impulse  which  will  result  in  a  given  deflection, When  more  exact 
experimental  values  are  available  for  the  applied  impulse  more  accurate 
computations  using  the  more  accurate  general  work  expression  of  Reference  1 
can  be  made. 

Bo  Buckling 

For  the  post  failure  buckling  region  the  circumferential  parameter  takes 
on  great  importance  and  the  simplification  given  by  equation  [12]  cannot  be 
employed.  Instead,  the  plastic  work  has  to  be  computed  from  the  more  gen¬ 
eral  integral  expressions  given  in  Reference  1.  We  use  equation  [40]  in 
the  "Errata  and  Addendum  of  Reference  1  and  assume  a  post  failure  deflec¬ 
tion  pattern  of  the  form 

c<4>*rr 

^  [23] 

=  ^v7 tjc  £  Ctrl'*. (irr~cf)  ir*<p<i.rr 

and  an  impulse  distribution  the  same  as  before,  i.e. 


-8- 


T  =  (ze)e~M  (/g*i) 


The  shell  used  was  made  of  aluminum  and  was 

in  diameter,  9"  long,  and  .006"  thick.  A  yield  stress  of  15,000  psi  was 
assumed.  Numerical  integration  of  equation  [40]  of  Reference  1  was  per¬ 
formed  assuming  a  perfectly  plastic  material  and  h>  ‘rX  =  5.  The 
value  of  '7L  =  5  was  obtained  from  equation  [10]  and  compared  well  with 
the  experimental  result.  The  calculated  peak  impulse  per  unit  area  to 
produce  ~  ,  S~  is  21  psi  milliseconds.  The  experiments  showed 

that  the  blast  load  which  produced  an  incident  impulse  of  21.1  psi  milli¬ 
seconds  and  reflected  impulse  of  46.8  gave  'ZO.Sm. 

VI .Conclusions 

One  cannot  apriori  assume  the  type  of  deformation  into  which  a  shell 
will  deform  under  blast  load  unless  some  pre-failure  computations  are  made 
to  determine  whether  the  she.1 1  will  collapse  or  buckle.  Schuman's  experi¬ 
ments^  point  out  that  these  two  types  of  patterns  can  exist.  The  compu¬ 
tations  and  experimental  comparisons  given  in  this  report  indicate  that 

1.  It  can  be  predicted  whether  a  shell  will  fail  in  buckling  or 
collapse  by  use  of  equation  [8]  and  [9]  of  this  report. 

2.  Once  the  failure  type  is  known  the  order  of  magnitude  of  the 
impulse  and  energy  to  produce  a  given  plastic  deformation  can  be 
computed  for  side-on  lateral  loads  by  assuming  a  collapse  pattern 
of  the  form  of  equation  [13]  for  collapse  failure  and  a  buckled 
pattern  of  the  form  [23]  for  post  buckling  behavior  in  the  plastic 
region. 

The  equations  given  in  Reference  1  have  been  programmed  for  the  BRL  com¬ 
puter.  This  will  enable  future  computations  with  more  complicated  de¬ 
flection  patterns. 
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